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We study the resonant inelastic x-ray scattering (RIXS) at the K edge of Ge. We measure 
RIXS spectra with systematically varying momenta in the final state. The spectra are a 
measure of exciting an electron-hole pair. We find a single-peak structure (except the elastic 
peak) as a function of photon energy, which is nearly independent of final-state momenta. We 
analyze the experimental data by means of the band structure calculation. The calculation 
reproduces well the experimental shape, clarifying the implication of the spectral shape. 
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1. Introduction 

Inelastic x-ray scattering is a promising method to study electronic structures in matters. 
Since intensities are much smaller than those on the elastic part, high-brilliance synchrotron 
sources are necessary. It is advantageous to use a resonant enhancement by tuning photon 
energy near absorption edge. This resonant inelastic x-ray scattering (RIXS) is described by 
a second-order optical process, in which a core electron is excited by an incident photon and 
then this excited state decays by emitting a photon to fill the core hole. Thereby electron- 
hole pairs remain in the final state. We can cover the momentum of the final state over the 
wide range by using the K edge of transition metals. Such RIXS experiments have recently 
been carried out on cuprates^'^ and manganites,^ having revealed characteristics of charge 
excitations through the momentum dependence of spectra. 

On semiconductors and insulators such as silicon,^ graphite,^ boron nitride,^ and others,^ 
RIXS experiments have already been performed, where corresponding photon energies are in 
the soft x-ray region. The spectra are a measure of exciting an electron-hole pair. Band struc- 
ture calculations have reproduced well the spectral shape due to small electron correlations. 
In this paper, we carry out a RIXS experiment using the K edge in Ge. Since corresponding 
photon energies are in the hard x-ray region, we expect that a new insight would come out 
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through the momentum dependence of the spectra. A similar experiment has aheady been 
reported.^ The present experiment is much thorougher with systematicaUy varying momenta 
in the wide range. Against our expectation, the experimental spectra as a function of photon 
energy show a broad single peak (except the elastic peak), and keep almost the same shape 
with varying momenta. 

We analyze the experimental data by calculating the spectra within the local density 
approximation (LDA). We disregard the Is core-hole potential in the intermediate state, 
expecting that this neglect would cause only a minor error. ^ This is partly justified from our 
calculation of the absorption coefficient, which agrees near the edge with the experimental 
one, as shown below. The final state consists of one electron in the conduction band and one 
hole in the valence band. If the core-hole level in the intermediate state is sharp enough, we 
can select the momentum of the excited electron in the conduction band by sharply tuning the 
incident-photon energy, because the denominator in the formula of the second order process 
(eq. (1)) selects a particular process through enhancement. Then, selecting the final-state 
momentum by setting the scattering geometry, we can specify the momentum of the hole in 
the valence band. In this situation, one expects sharp peaks at the energy of the electron-hole 
pair, as a function of photon energy. On the other hand, if the core-hole level is broad, many 
channels are opened to excite an electron with different momenta in the conduction band. We 
need to sum up each contribution. We calculate the RIXS spectra with varying values of the 
core-level width. This leads to an examination of determining the band structure from RIXS 
spectra. The l.s core level of Ge is rather broad with width of about 2 eV. Summing up 
each contribution of band-to-band transitions, we obtain the spectra in good agreement with 
the experiment. Thus the analysis based on the band structure calculation demonstrates the 
origin of the spectra. 

This paper is organized as follows. Experimental details are described in §2, and the for- 
malism for calculating the spectra is described in §3. In §4, the experimental data are presented 
and compared with the calculated results. Section 5 is devoted to concluding remarks. 

2. Experimental Measurements 

The experimental set up is shown in Fig. 1 (a). The geometrical arrangement {u!,9) is 
determined as shown in Fig. 1(b), for the momentum transfer from the incident X-ray with 
momentum to the outgoing X-ray with momentum q2. Experiment was performed at 
the bending magnet beam-line (7C) of the Photon Factory, Institute of Materials Structure 
Science. This beamline consists of a Si (111) double-crystal monochromator with the sagittal 
focusing mechanism and a focusing double mirror, and provides linear polarized x-rays. The 
incident beam size at the sample position was about 1.5 mm in horizontal and 1 mm in 
vertical; the photon flux was estimated to be ~ Ix 10^^ photons/s. The dimension of the 
Ge single crystal was 10 mm x 10 mm x 5 mm. Emitted x-rays were analyzed by means 
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of a cylindrically-bent InSb(880) crystal and measured by a position sensitive proportional 
counter. The total energy resolution of the emission spectrum in the present apparatus was 
3.0 eV. More details of the spectrometer has already been described in ref. 11. 

3. Formulation 



3.1 Double Differential Scattering Cross-Section 



RIXS is described by a second-order process. Using the generalized Fermi's golden rule 
where the interaction between light and electrons is treated by second order perturbation 
theory, we obtain the double differential scattering cross-section. 
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Here T = ei • Y^- ■pjC^'^^'^^ , T' = 62 • Pje~^^'^'^^ and e^, huji and hqi denote the polarization 
vector, the energy and the momentum vector of the incident (i = 1) and emitted (i = 2) pho- 
ton, respectively, and rj,Pj denote the position and momentum of electron j. Kets \i >, \m > 
and 1/ > represent initial, intermediate and finale states of the material system, respectively, 
with their energies Ei, Em, and Ej. The h/Tm is the life time of the intermediate state \m >. 

For germanium, an independent particle treatment for electron system seems to be ap- 
propriate, since electron correlations are weak.-'^'' For usual situations, the initial state \i > 
is assumed to be the ground state \Q >. The incident photon with momentum qi excites an 
electron from the Is core state c with crystal momentum kc to the conduction band (whose 
band index is e) with crystal momentum k(= qi + kc). Neglecting the core-hole potential 
in the intermediate state, we express the intermediate state by \m >= aj^gOk^^d^^ >• where 
a^(ag) is the crcation(annihilation) operator for the electron of ^ state. Then, the intermedi- 
ate state decays with emitting photon with momentum q2 and annihilating an electron in the 
valence band (whose band index is h) with crystal momentum k'(= q2 -|- kc) to fill the core 
hole. Therefore, we have the final state \ f >= aj^gak',/i|^ >• Since the final state uniquely 
determines the intermediate state, eq. (1) is rewritten as 
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where hijj = hui — huj2, and f2 = qi — q2 are energy and momentum of the final state. ee(k), 
e/i(k') and ec are the energy of the excited electron with crystal momentum k in the conduction 
band, that of the hole with crystal momentum k' in the valence band, and that of the core 
state, respectively. The crystal momentum conservation for the whole process is contained in 
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the factor of Kronecker 6, 
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where G is the set of reciprocal lattice vectors. Overlined quantities indicate their complex 
conjugates. Within the dipolc approximation, ta(p, r|e) = J dr^/;p^-r(r)e • p(t>a^{r — Ra), where 
R-o) V'pjT and (ji),^/''' arc the position vector of atom a in unit cell, Bloch-wavc function of a electron 
in the r-band with crystal momentum p, and Is-atomic orbital, respectively. These quantities 
are evaluated on the basis of the band structure calculation described in the following sub- 
section. 

3.2 Band Structure Calculation of Ge 

We perform a band structure calculation using the full-potential linearized augmentcd- 
planc-wavc (FLAPW) method"*^^ within the LDA. The local exchange-correlation functional 
of Vosko, Wilk and Nusair is employed. ""^^ The angular momentum in the spherical- wave ex- 
pansion is truncated at /max = 6 and 7 for the potential and wave function, respectively. The 
energy cutoff of the plane wave is 12 Ry for the wave function. It is well known that the 
density-functional theory underestimates bandgaps. For Ge, the density-fuctional calculation 
fails to obtain a finite bandgap: the valence-band top and the conduction-band bottom over- 
lap in energy at the F point, leading to a negative value of bandgap. Experiments suggests a 
finite bandgap, 0.744 eV. Many-body theories based on the Green function such as the GW 
calculation may be necessary to reproduce a collect size of the band gap. For the present 
purpose, however, it is sufficient to employ a simple remedy for the problem: carrying out a 
calculation on a system of a smaller lattice constant. The bandgap becomes larger for smaller 
lattice constants, since the splitting between bonding and antibonding bands becomes larger. 
We obtain a bandgap of 0.35 eV between the F and X points for the lattice constant a = 5.5A 
(a = 5. 658 A in the real system). Figure 2 shows the energy vs. momentum relation thus 
evaluated. The energy band is labeled by attached numbers for later use. Note that a slight 
change of the lattice constant has little influence on the dispersion relation except the band 



4. Comparison of experimental results and theoretical calculations 

We present the experimental data in comparison with the calculated spectra. The calcu- 
lation is made by following the procedure mentioned above with assuming F^ = F = 2.0eV 
(for all m) and Cc = —11, 103eV (measured from Fermi level). 

4-1 Absorption coefficient at the K edge 

Before going to the RIXS spectra, we first calculate the absorption coefficient at the K 
edge. Figure 3 shows the calculated result in comparison with experimental results by Kokubun 



gap. 
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et al}'^'^'^ The experiment shows no sign of core exciton. The calculated curve is in agreement 
with the experimental one in the energy region between 11,100 eV and 11,110 eV, indicating 
that the effect of core hole potential seems unimportant in the energy region of our interest. 
The calculation predicts extra peaks at about 11,115 eV and 11,125 eV, deviating from the 
experimental curve. This fault come from a linearizing procedure in the FLAPW method, 
which becomes less accurate for high energy states. 

4.2 RIXS spectra 

Figure 4 shows experimental RIXS spectra with varying incident photon energy between 
11,099 eV and 11,141 eV. The final-state momenta are set to be = —(3.5,3.5,3.5) and 
— (3.9,3.9,3.9). There are two peaks. One is an elastic peak, which is larger for the final-state 
momentum close to a reciprocal lattice vector. Another is an inelastic peak. Both get closer for 
the incident photon energy getting close to the threshold. Since the tail of elastic contribution 
masks the inelastic contribution, we have to subtract the former from experimental data in 
order to extract the latter contribution. Assuming that the elastic contribution dominates 
the intensity in the region oi fiuj = huji — hi02 < 0, we evaluate the clastic contribution for 
fiu) > Qhy the value at —huj. This may be a slight overestimate of the elastic contribution in 
the region of Juu > 0. Figure 5 shows the inelastic contribution evaluated this way in compar- 
ison with the calculation. The calculated intensity is normalized such that the peak height 
coincides with the experimental one at tuxJi = 11, 107 eV. We find a single broad peak, whose 
intensity increases with increasing incident photon energies. The calculation reproduces well 
such experimental behavior. This behavior of the peak intensity corresponds to the behavior 
of absorption coefficient at the K edge in the energy region between 11,100 eV and 11,107 eV 
where absorption coefficient increases with increasing incident photon energies (See Figure 3). 

Figure 6 shows RIXS spectra as a function of the Raman shift (= tkv) for fuvi = 11, 101 
eV. The final-state momenta vary from fi = —(3.1,3.1,3.1) to —(3.9,3.9,3.9). The elastic 
contribution is subtracted in the same way as above. The calculated results are in good agree- 
ment with the experimental ones. The spectra are nearly independent of final-state momenta. 
The reason will be discussed in the following subsection. 

The spectra are decomposed into each contribution of band-to-band transitions specified 
by band indices. Figure 7 shows the result for ft = —(3.3,3.3,3.3), Hcoi = 11,101 eV. We 
notice that the peak is mainly composed of transitions from the bands 3 and 4 to the band 5. 

4-3 Reconsideration of the idea of band- structure determination from RIXS spectra 

We examine more closely eq. (2) to clarify the origin of the spectral shape. This leads to 
the idea of band-structure determination from RIXS spectra proposed by Ma.^° 

We start by defining a function (here we neglect a-dependence of ta{p,T\e) and write 
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*(p,T|e)), 

^(k e) = u) 

(6e(k)-6,-/zu;i)2 + rV4- 

Since the core hole can carry momentum, the excited electron can take any value of momentum 
according to k = qi + kg. Here, we consider the following limit case opposite to the real 
situation in Ge. If the value of F is much smaller than the band-with of a branch e* in the 
conduction band, the denominator in eq. (4) can become very small for a certain value k* 
of k. In the case of band structure of Ge, the L point is selected when the incident photon 
energy is tuned at the absorption edge. Thus eq. (4) is approximated by 

V'(k,e) 0C(5(k-k*)(^e,e*- (5) 

Substitution of this relation into eq. (2) leads to 

oc ^ \tik* - n, h\e2)\' Siee* (k*) - 6,,(k* -n)-hw). (6) 

h 

The spectra consist of several (5-functions as a function of u, whose positions move with varying 
values of fl. Figure 8 illustrates how the band is determined when the incident photon energy 
is tuned at the absorption edge ( k* is at L-point). Thereby we can trace the valence band 
structure from the peak position by means of scanning ft. 

Next we consider the following limit case opposite to discussion mentioned above. For 
germanium, the value of F at the K-edge is estimated about 2 eV. This value is nearly equal 
to the band-with of the lowest conduction band of Ge. In such a situation, one may safely put 
the k-dependence on ■0(k, e), that is, put 

^(k,e) c^^(e). (7) 

This relation leads to an approximate expression to the differential cross-section, 

with 



^J2^ie)F{-re-hw), (8) 



F(^) = |«(k',^|e2)|'<5(^ - e^(kO), (9) 
(k',/») 

where is the average of ee(k) over k in the Brilloin zone. Since |i(k', /i|e2)P selects p 
symmetric states within the dipole approximation, F{Huj) is analogous to the density of states 
(DOS) projected onto p-symmetric states. The spectra are given by a superposition of the 
p-DOS, and this naturally explains why the spectra are nearly independent of final-state 

momenta. 

In Figure. 9, we show F-dependence of RIXS spectra. From above discussion, the spectra 
for F = 0.0001,0.01, (0.1) eV should be interpreted according to eq.( 6) and the spectrum for 
F = 2 eV should be interpreted according to eq.( 8). 
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5. Concluding Remarks 

Wc have carried out a RIXS experiment at tlie K edge of Ge witli systematically varying 
final-state momenta. We have obtained a broad inelastic peak as a function of photon energy, 
whose shape is nearly independent of final-state momenta. We have analyzed the experimen- 
tal data on the basis of the band structure calculation, summing up each contribution of 
band-to-band transitions. The band calculation has reproduced well the experimental data, 
demonstrating the origin of the spectra. 

We have examined a possibility of determining band structure from RIXS experiment. For 
smaller values of core-level widths, we can trace the energy dispersion with varying final-state 
momenta. To obtain sizable momentum change, however, photon energy is required to be 
larger than 2-3 keV. But, since corresponding core-level widths are usually large, the spectra 
are likely to be independent of final state momenta, thus making it difficult to determine the 
band structure. 
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Fig. 1. (a)Experimental set up for scattering, (b) Geometrical relations between qi,q2 and crystal 
axis of Ge. 
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Fig. 3. Absorption coefficient calculated at the K edge in comparison with the experiment (Ref.l5). 



11/ 



J. Phys. Soc. Jpn. 



Full Paper 



"(0 

c 



1 1 1 1 

(3.5 3.5 3.5) #1 


1 1 

\ 


11099 eV / 


\ 


11101 eV .....-•■^"^--A 
11103 eV 


1 ^B^'""»'^'**'t''»*V*»».«»„„.m>.f.«<n"n"-. 


mosev ^.--// 

11107 eV .--«»/^/ 




1 1 1 1 1 eV .^..x^ / 
11141 eV ^ 


» ^ A ^^^^^^^^^^ 


1 1 1 1 1 1 



11050 11060 11070 11080 11090 11100 11110 11120 
Emitted Photon Energy (eV) 



£1 

4—1 

"(0 

c 
a> 

c 



(3.9 3.9 3.9) 



' ■ tMK I ■ ■ • ■ • »W Ak« »'"4»%i*»«Hlt >1 1 



11099 eV 
"11101 eV„..^-^~0 



jmoT^ 



V" ■ W*'! ■•4( 



. -. ..11109 eM... 
^11111 gV 

11141 eV .....^^^^^^ 




11050 11060 11070 11080 11090 11100 11110 11120 

Emitted Plioton Energy (eV) 

Fig. 4. Experimental RIXS spectra as a function of emitted photon energy huji at (a) final-state 
momentum = —(3.5,3.5,3.5) and (b) —(3.9,3.9,3.9). The incident photon energy hjj\ varies 
from 11,099 eV to 11, 141 eV. 
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Fig. 5. Inelastic contribution of RIXS spectra as a function of emitted photon energy hu)2- Incident 
photon energy fkvi varies from 11,099 eV to 11, 107 eV. 



13/?? 



J. Phys. Soc. Jpn. 



Full Paper 





-20 -10 
-Raman Shift(eV) 



20 




(3.7 3.7 3.7) 



20 -10 
-Raman Shift(eV) 



10 



20 




(3.3 3.3 3.3) 



20 -10 
-Raman Shift(eV) 



10 



20 




-20 -10 
-Raman Shift(eV) 



20 



Fig. 6. Inelastic contribution of RXS spectra as a function of Raman shift fkj 
varying final-state momenta. 
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Fig. 7. Decomposition of RIXS spectra into contributions of band-to-band excitations, for tkJi = 
11, 101 eV, ri = —(3.3, 3.3, 3.3). Attached number m — n represents a transition from band m to 
band n. 
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Fig. 8. RIXS spectra calculated by assuming F = 0.0001 eV. 
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